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The Sagnac effect in two dimensional (2D) resonant microcavities is studied theoretically and 
numerically. The frequency shift due to the Sagnac effect occurs as a threshold phenomenon for the 
angular velocity in a rotating microcavity. Above the threshold, the eigenfunctions of a rotating 
microcavity become rotating waves while they are standing waves below the threshold. 
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The Sagnac effect is the phase difference between two 
counter-propagating laser beams in the same ring res- 
onator due to rotation, originally introduced by Sagnac 
in 1913[l|. It has become the basis for the operation of 
the optical gyroscopes such as ring laser gyroscopes and 
fiber optic gyroscopes after the invention of lasers and op- 
tical fibers in 1970's 2, 3, 4, 5] because the phase and fre- 
quency difference between clockwise (CW) and counter- 
clockwise (CCW) propagating beams are proportional to 
the applied angular velocity. These optical gyroscopes 
are normally used in airplanes, rockets, and ships etc. 
since they are the most precise rotation velocity sensors 
among any other types of gyroscopes. 

The Sagnac effect had been theoretically derived for 
the slender waveguides like optical fibers or the ring cav- 
ities composed of more than three mirrors by assum- 
ing that the light propagates one-dimensionally and the 
wavelength of the light is much shorter than the sizes 
of the waveguides or the rins; cavitiesHHI^. However, 
the sizes of the resonant cavities can be reduced to the 
order of the wavelength by modern semiconductor tech- 
nologies iSi? lE]- The conventional description of the 
Sagnac effect is not applicable to such small resonant 
microcavities. Especially, the resonance wave functions 
are standing waves which can never be represented by 
the superposition of counter-propagating waves. The as- 
sumption of the existence of CW and CCW waves plays 
the most important role for the conventional theory of 
the Sagnac effect. 

In this Letter, by perturbation theory typically used in 
quantum mechanics, we show that the Sagnac effect can 
also be observed even in resonant microcavities if the an- 
gular velocity of the cavity is larger than a certain thresh- 
old where the standing wave resonance function changes 
into the rotating wave. It is also shown that numerical 
results of the quadrupole cavity correspond very well to 
the theoretical prediction. Theoretical and numerical ap- 
proaches shown in this Letter do not assume that the CW 
and CCW waves exist in the cavity, but the pair of the 
counter-propagating waves is automatically produced by 
mixing the nearly degenerate resonance wave functions 
due to rotation of the cavity. 

According to the general theory of relativity, the elec- 
tromagnetic fields in a rotating resonant microcavity are 



subject to the Maxwell equations generalized to a non- 
inertial frame of reference in uniform rotation with angu- 
lar velocity vector 11 1, 2, 6, 10]. By neglecting 0(/i^), we 
obtain the following wave equation for the electric field 
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where h = ^ {r x ft) . In the above, c and n are respec- 
tively the velocity of light and the refractive index inside 
the cavity. 

In conventional theoretical approach for the Sagnac ef- 
fect, the frequency shift of the resonance proportional to 
the angular velocity of the rotating ring cavity is de- 
rived from assuming that the electric field in Eq. (1) 
propagates one-dimensionally along the slender optical 
waveguides 0, 0. This method is not applicable to the 
resonant microcavities because the wavelength of the res- 
onance is not much shorter than the size of the cavity and 
the electric field does not propagate one-dimensionally. 

Instead, in the case of the 2D resonant microcavity per- 
pendicular to angular velocity vector fi, the resonances 
can be obtained by solving the following stationary wave 
equation derived from Eq. (1) for the stable oscillation 
solution. 
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where the 2D resonant cavity is rotating on xy-pldine 
clockwisely,i.e., ft = (0,0,^1) and Q > 0. We as- 
sumed that TM mode of the electric field oscillates as 
E{r, t) = (0, 0, V^(r)e-^^^^ + c.c). 

For simplicity, we impose the Dirichlet boundary con- 
dition on the electric field of the resonant microcavity in 
the remainder of this Letter. 

In the case of a microdisk cavity, Eq. (2) can be solved 
exactly as follows. Eq. (2) is rewritten in the following 
form in the cylindrical coordinates. 
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One can assume the solution ?/;(r, 6) is given as ^|){r^ 6) = 
/(r)e*^^ where m is an integer, and then obtains 
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f{r) = 0, 
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Eq. is the Bessel differential equation, and so the solu- 
tion /(r) should be the Bessel function of the m^^ order 
Jm{Kmr). The eigenvalue of the wave number k is given 
by the zero of Jm{KmK) because of the Dirichlet bound- 
ary condition, where R is the radius of the microdisk 
cavity. 

Accordingly, the shifted wave number due to rotation 
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where ko is the zero of the Bessel function Jm{nkR), that 
is, the eigenvalue of the wave number when the cavity 
is not rotating. Consequently, when the microdisk cav- 
ity is rotating, the wave function is the rotating wave 
Jm{Kmr)e^'^^ and the degenerate wave number k^ with- 
out rotation splits into two different wave numbers of the 
counter-propagating waves corresponding to the signs of 
the integer m. From Eq. 0, one obtains the frequency 
difference Acj between the counter-propagating waves. 
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It is important that the CW and CCW wave solutions 
are degenerate eigenstates even when the microdisk cav- 
ity is not rotating, and hence the standing wave solution 
produced by the superposition of these degenerate rotat- 
ing waves are also the eigenfunction. However, with a 
finite angular velocity 1^, the CW and CCW wave solu- 
tions become non-degenerate states, which means that 
only the rotating waves are the eigenfunction of the ro- 
tating microdisk cavity. The frequency difference Ao; 
between the CW and CCW solutions is proportional to 
the angular velocity ^. This is the Sagnac effect for mi- 
crodisk cavities. 

In general cases of 2D resonant microcavities of arbi- 
trary shapes, Eq. (2) cannot be directly solved in the 
same way for microdisks. 

First, we discuss the case that the spacing A/c between 
the adjacent eigenvalues k of the wave number is large 
enough to satisfy the following inequality. 
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where ipo and ipi are the wave functions of these eigen- 
states which correspond to the adjacent eigenvalues when 



the angular velocity Q is zero. We assume that, due 
to the rotation of the cavity, the eigenvalue is shifted 
as /c = ko -\- 5k and the wave function is changed as 
= V^o + Y.1^0^1^1^ where (V^ ^n^ki'^) ipi = 0. Here 
Sk and ci are assumed to be so small as Q/c. Then, from 
Eq. ^ we obtain 
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Using the following relation, 



V^oV 'ds = ^, (10) 



where D and dD denote respectively the domain and the 
edge of the cavity, we finally obtain ^/c = up to the first 
order of |r^/c|, and 
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Consequently, as long as the angular velocity is small, 
there is no Sagnac effect, which means that the 
wave functions are standing waves instead of counter- 
propagating waves and the frequency difference between 
two standing waves does not increase. 

Next we discuss the case that the spacing Ak between 
two wave numbers ko and ki is so small that it does 
not satisfy the inequality According to the pertur- 
bation theory for nearly-degenerate states in quantum 
mechanics, the wave function should be represented as 
the superposition of two nearly-degenerate eigenfuncions: 
ip = coV^o + cii/ji. Substituting this equation into Eq. ((^J 
yields 
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where M is the following matrix: 



n^{P - ko^) -2ik J drtlJo (h • V) V^i 

-2ik J drilJi {h • V) V^o n^{k'^ - h^) 

In order to obtain non-trivial solutions for Eq. ([T^ . the 
determinant of M should vanish, which yields a quadratic 
equation for k'^. Consequently, we obtain the eigenvalues 
of the wave number up to the first order of |^/c|. 
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Accordingly, the frequency difference Auj between the 
two eigenfunctions newly produced by rotation of the 
cavity is proportional to the angular velocity llj. 



Acj 
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Then, from Eq. (|T^ we also obtain the ratio of the coef- 
ficients Co and ci as follows: 



(15) 



It is important that the standing wave function tpQ 
without rotation of the cavity can be represented by 
the superposition of the Bessel functions as i/jq = 
"^"^^Q am Jmi'^kor) cos mO , if the cavity is symmetric 
for X-axis and tpo is an even function. The other wave 
function i/ji nearly degenerate to i/jq which should be 
an odd function for x-axis also can be written as tpi = 
Z)m=o bmJm{nkir) sin mO , where ^ bm and ko ^ ki. 
When the cavity is rotating but the angular velocity is 
small enough to satisfy the inequality Q, the frequency 
difference does not increase and the wave functions do 
not change drastically. However, when the angular ve- 
locity is increased to be larger than this threshold which 
violates the inequality ((Hj), the wave functions change into 
the CW and CCW rotating waves because of Eq. (|T5|) . 
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m=0 

CO 

= J2''rnJm{nkr)e^'"'', (16) 

m=0 

where the wave number k is given by Eq. (O. Then, 
according to Eq. (|T^ . one can observe the Sagnac ef- 
fect as the frequency difference Auj between counter- 
propagating waves proportional to the angular velocity. 
The most important point is that there is a threshold 
of the angular velocity for the Sagnac effect in resonant 
microcavities. It is possible to delete this threshold by 
symmetric shapes of the cavities as shown in the case of 
microdisks. 

Next we numerically solve the Maxwell equation for a 
rotating microcavity, and show that the threshold phe- 
nomena of the transition from the standing wave solu- 
tions to the counter-propagating wave ones actually oc- 
cur. 

First, let us explain the numerical method briefly. The 
solutions '0(r, 0) for the Maxwell Equation for the 
stationary solutions of the rotating cavity can be repre- 
sented by the superposition of the Bessel functions in the 
cylindrical coordinates. 
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Then the Dirichlet boundary condition is imposed as 
iIj{R{6)^6) = 0, where R{6) denotes the edge of the cav- 
ity. Accordingly, we obtain 
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For numerical computation, the infinite sum over m can 
be approximated by the finite sum from — M to M where 
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FIG. 1: The nearly- degenerate wave functions of the 
non-rotating quadrupole cavity corresponding to the eigen- 
wave numbers {a)nkARo = 49.3380585 and (h)nkBRo = 
49.3380615. We call the modes of (a) and (b) modes A and 
B, respectively. 



M is a large integer. Then Eq. ([T^ can be rewritten as 
the matrix form. 
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0, (19) 



where Anm = J^"" Jm{KmR{0))e'^'^-''^^dO. Therefore, 
the eigenvalues of the wave numbers can be obtained nu- 
merically as the zeros of the determinant of the above 
matrix because of non-trivial solutions for the coefficients 

For numerical calculation, we choose a quadrupole cav- 
ity, which is defined by the boundary R{0) = Ro{l -\- 
ecos26>)0. The parameters of the quadrupole are set 
as follows: Rq = 6.2866/im, e = 0.12, and the refractive 
index n = 1. When the quadrupole cavity is not rotat- 
ing, solving the Helmholtz equation Eq. (3) with 1^ = 
yields the nearly-degenerate standing wave eigenfuctions 
as shown in Fig. 1. We call the two modes shown in 
Fig. 1(a) and (b) modes A and B, respectively. 

When the angular velocity Q is smaller than a certain 
threshold Qth (where RoQth/c = 5.0 x 10~^), the fre- 
quency difference does not increase as shown in Fig. 2(a), 
and the eigenfunctions remain standing waves. However, 
for ft > Vtthi the frequency difference increases gradu- 
ally, and becomes proportional to Vt^ and modes A and 
B drastically change into the rotating wave functions as 
shown in Fig. 3(a) and (b), respectively. 

The transition from the standing wave to the rotating 
wave can be clearly observed also by the CW and CCW 
wave components icw.ccw defined as follows based on 
rotating wave decomposition of Eq. (|T7|) : 



icw{ccw) 



(20) 



m<0(m>0) 



where ^ Zlm/o l^^nP- For < Vtth, icw.ccw are 
around 0.5 for modes A and B as shown in Fig. 2(b) and 
(c). When Vt exceeds ^th^ iccw{cw) of mode A(B) sud- 
denly vanishes, which means the wave function consists 
of only CW(CCW) waves. Therefore, one can see that 
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FIG. 2: (a)The (dimensionless) frequency difference RqAco/c 
versus the (dimensionless) angular velocity. The frequency 
difference does not change for RqQ/c < RoQth/c{= 5.0 x 
10~^). For RqQ/c > RoQth/c, the frequency difference be- 
comes proportional to the angular velocity. The CW and 
CCW waves components versus the angular velocity are 
shown for (b) mode A and (c) mode B. 

modes A and B in Fig. 3 are the CW and CCW rotating 
waves, respectively. 

In summary, we have shown that the Sagnac effect can 
be observed in 2D resonant microcavities when the angu- 



[1] E. J. Post, Rev. Mod. Phys, 39, 475 (1967). 

[2] W. W. Chow, J. Gea-Banacloche, L. M. Pedrotti, V. E. 

Sanders, W. Schleich, and M. O. Scully, Rev. Mod. Phys. 

57, 61, (1985). 

[3] F. Aronowitz, in Laser Applications, M. Ross, ed. (Aca- 
demic, New York, 1971), Vol. 1, pp 133-200. 

[4] Vah V. and R. W. Shorthill, Appl. Opt., 15, 1099, (1976). 

[5] S. Ezekiel and S. R. Balsamo, Appl. Phys. Lett. 30, 478, 
(1977). 

[6] L. N. Menegozzi and W. E. Lamb, Jr. Phys. Rev. A, 8, 
2103, (1971). 

[7] Y. Yamamoto and R. E. Slusher, Physics Today 46, 66 
(1993). 

[8] R. K. Chang and A. J. Campillo, eds.: Optical processes 



lar velocity is larger than a certain threshold where the 
nearly degenerate standing wave eigenfunctions of the 
non-rotating cavity change into the pair of the counter- 
propagating waves. 

The threshold phenomenon seems to be akin to the 
lock-in phenomenon, which occurs owing to the mode- 
locking between the counter-propagating waves 2, "s*, 
However, in our theoretical approach, the effects of 
backscattering and an active medium which both cause 
the lock-in phenomenon are not taken into account. The 
existence of the threshold can be shown even without 
these effects. 

The Sagnac effect in microcavities will be actually ob- 
servable by measuring the frequency difference in a ro- 
tating frame of reference, as the measuring method of 
optical gyroscopes 2, 3]. A discussion on the actual ex- 
periment will be reported elsewhere. 




(a) (b) 

FIG. 3: The wave functions of the rotating quadrupole cavity 
with the (dimensionless) angular velocity RqQ/c = 6.28 x 
10~^(> Ro^th/c) respectively corresponding to (a)mode A 
and (b)mode B. 
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